Ferromagnet-superconductor proximity effect: The clean limit 
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We study theoretically the influence of ferromagnetic metals on a superconducting film in the 
clean limit. Using a self-consistent solution of the Bogoliubov-de Gennes equation for a ferromagnet- 
superconductor-ferromagnet double junction we calculate the pair potential and conductance spectra 
as a function of the superconducting layer thickness d for different strengths of ferromagnets and 
interface transparencies. We find that the pair potential and the critical temperature are weakly 
perturbed by the exchange interaction and do not drop to zero for any finite d. On the other 
hand, for thin superconducting films charge transport is spin polarized and exhibits a significant 
dependence on the ferromagnetic strength and magnetization alignment. 



I. INTRODUCTION 

Rapid advancement of nanofabrication technology in 
the past decade has reinvigorated interest in under- 
standing the effects inherent to clean superconducting 
heterostructures. Apart from potential device appHca- 
tions, particularly in quantum information storage and 
processing/ variety of phenomena makes systems con- 
sisting of superconductor sandwiched between two ferro- 
magnets interesting for both experimental and theoret- 
ical investigation (see, for example, i-efsi2i2iiiSiSiL& and 
9,10,11,12,13,14,15,16,17,18,19,20)^ The interplay of ferromag- 

netism and superconductivity results in characteristic 
proximity effects near the contacts between two metals: 
the Andreev reflection is suppressed in the presence of 
a ferromagnet^ conversely, the superconducting corre- 
lations, described via the pair amplitude, extend in the 
ferromagnetic material in an oscillatory way^ 

The proximity effect in heterostructures consisting of 
a superconductor (S) in contact with a normal, nonmag- 
netic (N) or ferromagnetic (F), metal has been widely 
studied and well understood in the dirty limit, where the 
electron mean free path I is much smaller than the su- 
perconducting coherence length In particular, 
a strong depairing effect is found in FSF trilayers, which 
results in a superconducting-to-normal phase transition 
at a finite thickness of the S layer Furthermore, using a 
spin-valve setup for these junctions, spontaneous transi- 
tion from parallel (P) to antiparallel (AP) ferromagnetic 
moment orientation could be triggeredri^ 

However, in recent experiments with FSF trilayers it 
was found that the difference between the superconduct- 
ing critical temperatures for AP and P alignment is less 
by two orders of magnitude than theoretical predictions 
for the dirty limit >^ A possible reason for this disagree- 
ment could be that these samples were actually much 
cleaner than the ones used in previous experimentsiSi^ 
Namely, as a consequence of phase coherence of electron 
wave functions quasiparticle spectrum of clean trilayers 
differs substantially from the BCS result for bulk super- 
conductor: density of states is practically gapless and 
Andreev reflection is reduced in thin S layers. ^^'^^ As a 
result, depairing induced by ferromagnets appears to be 
much weaker in the clean limit, when Z ^ 



The purpose of this paper is to clarify the ferromagnet- 
superconductor inverse proximity effect in the clean limit, 
within the framework of the BCS theory. The model 
we study is an FISIF heterostructure with an S film in 
contact with massive F metals, where I denotes the in- 
terface potential barrier of arbitrary transparency. To 
describe the equilibrium and transport properties of 
such a system, we use a numerical procedure to solve 
the Bogoliubov-de Gennes equation self-consistcntly for 
spatial-averaged pair potential. We find that neither the 
magnitude of exchange interaction in the ferromagnets 
nor the relative orientation of the magnetizations has any 
significant infiuence on the pair potential and the criti- 
cal temperature of the superconductor. On the other 
hand, transport properties of thin S films exhibit a sig- 
nificant dependence on ferromagnetic strength and align- 
ment of magnetizations. ^^•^'^ Unlike the case of diffusive 
mesoscopic superconducting bilayers and trilayers jl5 we 
found a BCS type of the pair potential temperature de- 
pendence. A peculiar temperature and alignment depen- 
dence of the proximity effect is recently obtained for a 
model of FSF trilayers with atomic thickness iSSiSi 



II. EQUILIBRIUM PROPERTIES 

In this section we study the equilibrium properties of 
an s-wave superconductor sandwiched in between two fer- 
romagnets. We consider a simple model of FISIF double 
junction consisting of a superconducting layer of thick- 
ness d connected to ferromagnetic metals by interfaces of 
arbitrary transparency, fig. ^ Assuming that the metals 
are clean, quasiparticle propagation is described by the 
Bogoliubov-de Gennes equation 
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where i?o(r) = -h^^"^ /2m + W{v)- Ep, Ep is the Fermi 
energy, cr is the quasiparticle spin orientation (tr =t,i 
and a and E is the quasiparticle energy with 

respect to the Fermi level. The interface potential is 
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modeled by W(r) = W{S{z) + S{z — d)}, where the 
z axis is perpendicular to the layers and S{z) is the 
Dirac delta- function. For simplicity, the ele ctron effec- 
tive mass m and the Fermi velocity, vp = \/2mEp/h'^, 
are assumed to be constant through the junction. The 
superconducting pair potential is taken in the form 
A(r) = A(z)8(z)8(d — z), where 0(z) is the Heaviside 
step function. The exchange potential h{r) is given by 
ho{Q{~z) + [~]<d{z - d)} for the P [AP] alignment, and 
Per is 1 (—1) for a =t (i); a uniform magnetization is 
assumed to be parallel to the layers. The parallel com- 
ponent of the wave vector, q|| = {q^, 0), is conserved, 

and the spinor (M(j(r) Wg.(r))"^ satisfies the appropriate 
boundary conditions for the wave function and its first 
derivative at 2; = and z = di^ 



summation over q, we get 




FIG. 1: The geometry of FISIF heterostructure in P (AP) 
alignment of magnetizations. 

For a constant pair potential, solution of eq. ^ has 
the following form in the superconductor, < z < d, 



= cxp(jq|| • r) X 



(2) 



x| [ci{E,q^i)cxp{iq+z) + C2{E,q^i)cxp(-iq+z)] 
+ [c3{E,q\i)exp{iq-z) + Ci{E,qi\)eiip{-iq^z)] f ^ j } 



Here, u ^ ^(1 + n/E)/2 and v = ^(1 - n/E)/2 are the 
BCS coherence factors, ft = Vi?^ — A^, and the mod- 
ulus of the wave vector q^ = {qx,qy,qf) is given by 



■\/ {2m/h?){EF ± ri). Coefficients ci through C4 
are obtained from boundary conditions in the scattering 
problem for FISIF heterostructure.^'^ 

Wave-vector components have to be normalized by us- 
ing the condition for canonical transformation 



/ d'r[\u^ir)f + \v^{r)f] = 1, 
Jv 



(3) 



where the integration is performed over the volume V of 
the superconductor. The self-consistency condition for 
the pair potential is given by^S 

A(z) = Aj] {uT(r>*(r) [1 - f,{E)] - u^{r)v^r)fo{E)} . 
q 

(4) 

Here, A is the coupling constant and /o is the Fermi dis- 
tribution function at temperature T. By performing the 



A{z)^XN{0)V d6isin6'cos6' X 
"'0 



dn {u^{r)vl{r) [1 - fo{E)] - u^{r)v*{r)fo{E)} 



where -/V(O) = mkF/2TT'^h? is the normal-metal density of 
states (per spin orientation) at the Fermi level, 9 is the 
angle between q"*" and the z axis, and Hujd is the upper 
cutoff in integration over quasiparticle kinetic energy Q. 




FIG. 2: Zero-temperature pair potential A, normalized to 
the bulk value Aq, as a function of the S film thickness d 
for NSN (Z = 0) and NISIN (Z = 10) junctions. The same 
curves describe dependence of Tc/Tco on d. 

We solve eq. using the stepwise approximation for 
A(z). We calculate the spatial average of A(z) following 
a standard iteration procedure 



1 r 

A,+i = -J A,{z)dz 



(5) 



by setting E = a/TP + AI in eq. Q in order to obtain 
Ai(z) in the i-th iterationiS Starting from the bulk value 
Ao we repeat this procedure until the difference between 
Ai+i and Ai becomes sufficiently small. This procedure 
is justified both for thin S films, d/^o ^ I7 where the 
spatial variation of the pair potential is small, and thick 
S films, d/^o ^ 1, where A(z) is practically flat at Aq, 
except in a narrow region of the order of ^0 near the 
interfaces. 

Numerical results obtained by the self-consistency al- 
gorithm described above are shown in figs. [21 and |31 
The average pair potential A, normalized to the bulk 
value Aq, is shown in fig. [21 as a function of cJ/^Qj where 
^0 = fivp /t^Aq, for an NISIN junction at zero temper- 
ature and for two strengths of the interface barriers, 
Z = 2W/hvp. As the transmissivity of the interfaces 
decreases (i.e., as Z increases), the normal reflection of 
incoming electrons becomes more probable, which makes 
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FIG. 3: Zero-temperature pair potential A, normalized to the 
bulk value Aq, for (a) FSF {Z = 0) and (b) FISIF {Z = 10) 
junctions as a function of the S film thickness d for X — 0.5 
and for parallel and antiparallel magnetization alignments. 
The corresponding NISIN curves {X = 0) are shown for com- 
parison (solid lines). The same curves describe dependence of 
Tc/Tco on d. 



the S film more isolated and weakens the proximity elFect. 
This could be seen for Z = 10, which corresponds to a 
single-barrier transmissivity of ~ 1% of the electrons in- 
jected into the superconductor. In the case of transparent 
interfaces, Z = 0, the pair potential is more suppressed 
by the proximity of the normal-metal electrodes. 

We have studied ferromagnetic influence for various 
strengths of the normalized exchange potential X = 
ho/Ep in the F metals for both P and AP magnetiza- 
tion alignment In fig. O we show the particular case 
of X = 0.5, for FSF {Z = 0) and FISIF {Z = 10) het- 
erostructures. It can be seen that all the points, cal- 
culated by numerical integration, fall very close to the 
X = curves, almost independently of the magnetiza- 
tion alignment. This is in a strong contrast with dirty 
FSF structures, which exhibit a phase transition to nor- 
mal state at some critical value of the superconducting 
film thickness2i2. which varies with relative orientation of 
the magnetizations We conclude that in clean FISIF 
systems the superconducting pair potential is practically 
independent of ferromagnetic polarization. This result 
is in accord with the previous ones obtained by a fully 
self-consistent description of the proximity effect at the 



F-S interfaceJfi Recent experiments also imply that the 
difference between the superconducting critical tempera- 
tures in AP and P alignment is off by about two orders of 
magnitude with respect to the theoretical predictions for 
the dirty limit. ^-^ A possible reason for this disagreement 
could be that these samples were much cleaner than what 
was the case in earlier experiments (refsi^ and^) where 
data fit very well to the dirty-limit theory. 

Qualitatively, the fact that A weakly depends on X 
in the clean limit can easily be shown using the quasi- 
classical approximationi^iiS^ By neglecting the wave vec- 
tors outside the small interval around the Fermi surface 
of radius kp, so that h^q^/2m ~ —Ep ± ti^kplql/m, 
the one-electron Hamiltonian could be linearized, Hq ~ 
zLi{h'^kp/m)(d/dz) — 2Ep. After normalization given by 
eq. coefficients in eq. become 



Ci(S,q||) = 
C3(-E,q||) = - 



A/y(l-u2w2C-isinC)' 



C2(£;,q||) = 0, 
Ci{E,ci\\) = 0, 



where C, = d{qf — qj). Therefore, solutions ^^(r) and 
Va (r) of linearized eq. (Q) are independent of X and ex- 
actly match the solutions of the Bogoliubov-de Gennes 
equation for transparent NSN junctions. We emphasize 
that a weak influence of the exchange potential is inher- 
ent for clean S and F metals. The dependence of the criti- 
cal temperatures on the mutual orientations of ferromag- 
netic moments is also hardly observable in dirty FISIF 
hybrids with finite interface transparency^^ii^ However, 
for such structures with higher transparency the inverse 
proximity effect is significant. ^-^ 

We have also verified that the temperature dependence 
of A remains to be of the BCS type,^^ well described by 

A(T) = A(0)tanh (^l.Ti^/Tc/T - 1^ for aU d, where Tc 

is the critical temperature of the S film. In addition, the 
d-dependence of the critical temperature normalized to 
the bulk value, Tc/Tco, coincides with A(0)/Ao(0) vs. d 
curves (figs. [21 and O), being greater for AP than for P 
alignment of magnetizations, which is in agreement with 
theoretical predictions ^^-^^ and experimental results^.. In 
diffusive mesoscopic superconducting bilayers and trilay- 
ers, however, a more complex dependence is predictedi^ 
Moreover, for finite F layers the critical temperature, 
as well as the pair potential, oscillate with the F-layer 
thickness ,SiSiii 



III. TRANSPORT PROPERTIES 

The properties of quasiparticle transport are com- 
monly described by differential conductance spectrum. 
For an FISIF double-barrier junction at zero tempera- 
ture conductance per orbital transverse channel can be 



4 



1.0 



0.5 



0.5 



O.Q 



' 1 ' 


— p 


x=o 




-AP 

z=o 


X=Q.5 

i_l 1 1 1 1 ■ ■ ■ ■ 


Z=l 

— p 


1 ' 

^.I^----. A=0 


-AP 


X=0.5 



FIG. 4: Zero bias conductances of FSF {Z = 0, top panel) 
and FISIF {Z = 1, bottom panel) junctions at zero tempera- 
ture as a function of the S film thickness d for X = 0.5. The 
corresponding NSN and NISIN spectra {X = 0, dashed lines) 
are shown for comparison. 

calculated fromMJ^ 
G{E) =GoY.P-^ 

X / deisinffcosff [A„{E,9) +C„{E,9)]i6) 
Jo 

where Gq — 2e^//i is the conductance quantum and 
Pct = (1 + PaX) /2. Analytical results for the probabili- 
ties of Andreev reflection and direct quasiparticle trans- 
mission, Acr{E,6) and GaiE,9), respectively, are pre- 
sented and discussed in refii^. Contribution of evanes- 
cent propagation to G{E) is included as well. Here, we 
only point out the strong influence of ferromagnetism on 
charge transport through a thin S fllm, in contrast to the 
weak proximity effect on equilibrium properties. 

As an illustration, we calculate the zero bias conduc- 
tance, G(0), using the self-consistent pair potential ob- 
tained in the previous section. The results are shown 
in flg. ^ The subgap transport through the supercon- 
ductor changes by virtue of two principal mechanisms: 
firstly, due to the presence of the ferromagnets, and sec- 
ondly, due to the decrease of thickness d. The zero-bias 
conductance for an FSF junction with X = 0.5 is sig- 
nificantly below the NSN value of unit conductance per 
channel, and splits for the P and the AP magnetization 



alignment as the S film becomes thinner. The subgap 
conductance is greater in P than in the AP alignment 
as a result of a strong magnetoresistive effect in thin S 
layers: when d ^ the direct transmission of spin polar- 
ized quasiparticles across the superconductor becomes a 
dominant transport mechanism>iii^ Hence, in contrast 
to the NSN trilayers, the conductance of an FSF sys- 
tem shows a strong size effect. Conductances become 
significantly suppressed for finite interface transparency 
(Z = 1), both for NISIN and FISIF junctions, due to 
the increase of normal and decrease of Andreev reflec- 
tion probability. 



IV. CONCLUSION 

We have studied the ferromagnet-superconductor 
proximity effect in clean FISIF heterostructures with thin 
S layers and massive F metals. We have found that the 
superconducting order parameter is weakly affected by 
the exchange interaction and has practically the same 
dependence on the S fllm thickness as in the correspond- 
ing NISIN structures. On the other hand, quasiparticle 
dynamics within an FISIF heterojunction is substantially 
changed with respect to the corresponding NISIN system, 
the more so as the S fllm becomes thinner. 

While in clean ferromagnet-superconductor hybrids 
with massive F layers the ferromagnetism has negligible 
inverse proximity effect on equilibrium properties, the 
opposite behavior is previously found for the dirty ones;^ 
where the pairing potential may be very sensitive to the 
vicinity of the F layer(s). Moreover, in dirty hybrids the 
critical value of superconducting layer thickness at which 
destruction of superconductivity occurs is strongly de- 
pendent on the ferromagnetic exchange potential. On 
the other hand, charge transport is diffusive in this case, 
and consequently the conductance spectrum is practi- 
cally independent on the S film thickness and alignment 
of magnetizations. 

In summary, we have found that in the clean limit 
the BCS self-consistent solution for a thin supercon- 
ducting film in contact with massive ferromagnets shows 
very weak depairing effect of the exchange interaction. 
While the equilibrium properties are practically unaf- 
fected, charge transport is strongly influenced by prox- 
imity of ferromagnets. 
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